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Abstract 

We have studied numerically the frequency distribution p{n) of the n-th neighbour along the 
optimal tour in the Euclidean travelling salesman problem for N cities, in dimensions d = 2 
and d = 3. We find there is no significant dependence of p(n) on either the number of cities 
N or the dimension d. 
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1. Introduction 

The study of optimization problems is of considerable interest to computer scientists, math- 
ematicians and physicists alike. In a typical optimization problem, there is a large finite 
set of possibilities to search from, in order to obtain the optimal solution: if the problem 
is of "size" A^, then typically, there are of the order of A^! or possibilities, of which we 
want the one that minimizes (or maximizes) the cost function. The travelling salesman 
problem (TSP) is a simple example of a combinatorial optimization problem where, given 
a certain set of cities and the intercity distance metric, a travelling salesman must find the 
shortest tour in which it visits all the cities and returns to its starting point |[l|, H|. It is 
a non-deterministic polynomial complete (NP-complete) problem. There are two forms of 
TSP which are of interest: the Euclidean TSP and the random link TSP. In the Euclidean 
TSP, the N cities are distributed with uniform randomness in a (i-dimensional hypercube 
and the distance is measured in the Euclidean metric (with AZ = v^Ax^ + Ay^), whereas 
in the random link TSP, the distances between the cities i and j are taken as independent 
random variables with a given distribution. It was noted that the random link TSP can be 
mapped onto the Euclidean model, provided the distribution is chosen appropriately and the 
correlations between three or more distances neglected 0]. 

A city is said to be the n-th neighbour of a reference city if there are exactly (^ — 1) other 
cities that are nearer to the reference city. In a given configuration of cities, for every city 
we can find its neighbours, arrange them in order of their distances from that city and label 
them consecutively with their neighbour number n. Thus n = 1 is the nearest neighbour, etc. 
We can find out how many times the n-th neighbour is chosen along the optimal tour and 
determine the frequency distribution p{n). Here, we have studied numerically the frequency 
distribution p{n) along the optimal tour for the Euclidean TSP only. The optimum tours are 
obtained with the help of branch and bound algorithms with open boundary conditions. As 
a matter of interest, we have also studied the frequency distribution p(n) using the "greedy 
algorithm". 
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2. Simulation and numerical determination of p{n) 

We generated random configurations in dimensions d = 2 and 3, for diff'erent sizes = 10 
to 256. For each configuration, we determine the frequency distribution of neighbours p{n) 
along the optimal tour and then take the average over 5000 configurations. We have found 
the optimum tours m d = 2 and 3, with the help of branch and bound algorithms using open 
boundary conditions (Fig. 1). The results for p{n) in dimensions d = 2 and 3 are shown 
in Fig. 2 (a) and (b). The numerical values of p(n) for different values of n in d = 2, are 
plotted against in Fig. 2 (c). These show that the frequency distribution p(n) does not 
vary significantly with and so p{n) does not have any prominent finite-A^ effect, and it 
does not depend on the dimension d. The frequency distribution of neighbours, obtained by 
numerical fitting, is of the form p{n) = Aexp{—an)[l + B exp{n/b)] where A = 0.72 ± 0.01, 
a = 0.60 ± 0.03, B = 0.05 ± 0.01 and 6 = 3.0 ± 0.1. The errors in A, a, B and b are obtained 
by eye-estimation. 

As a matter of interest, we have also studied p(n) using the "greedy algorithm", where 
from each city the salesman goes to the nearest city not already in the tour and finally from 
the A^-th city returns directly to the first. The results for p{n) in dimensions d = 2 and 3 are 
shown in Fig. 3. Here also, p(n) does not depend on the dimension d. There seems to be a 
cross-over from the exponential decay to a power law decay at n = 6, for some unapparent 
reasons. 

3. Discussions and summary 

It was shown in ||l| that the average optimal travel distance in the unit d-dimensional 
hypercube, where A^ random cities are distributed uniformly, < /j^^ > scales as N^~^/'^. 
The empirically obtained frequency distribution p{n) along the optimal tour can be used 
to get a rough estimate of this average optimal travel distance < l^^^ >. We may write 

< zjy-* >~ Yln=i < P(^) X > where < ... > denote the ensemble averages and 

< D^j^\n) > is the average n-th neighbour distance along the optimal tour between N cities 



2 



in a unit d-dimensional hypercube. Note that the average optimal travel distance < /j^'' > 
should actually involve the average over the product < p{n)D^^\n) >, while we use here the 
product of averages < p{n) >< D^]^\n) >. Here, we use the general expression given in 
for the n-th neighbour distance in the unit d-dimensional hypercube containing random 
points distributed uniformly: < D'^^\n) >= [T{d/2 + 1)/ (n'^/^N)]^^'^ [T{n + l/d)/T{n)], and 
the empirically determined frequency distribution p{n) (=< p{n) >) to estimate < 1^^ >■ 
We find that for = 2, < 1^ >~ 0.77N^/^ and for d = 3, < 1^ >~ 0.74Ar2/3_ ^^^y 
noted that the ensemble of D^^\n) in TSP is not the same as the case where distances are 
calculated without any restrictions. 

In summary, here we study numerically the frequency distribution p{n) along the optimal 
tour for the Euclidean TSP for dimensions d = 2 and 3. For optimum tours, we find there is 
no significant dependence of p{n) on either the number of cities or the dimension d, and 
the empirically determined frequency distribution p{n) ~ Aexp{—an)[l + B exp{n/b)]. Since 
p{n) does not have any prominent finite-iV effect and the values of p{n) remain significant 
for small n (up to n ~ 45 in our study for < 256), one can therefore determine p(n) quite 
accurately using small system sizes, and thus optimizing the computational efforts. 

Acknowledgement : The author is indebted to Prof. B. K. Chakrabarti for introducing 
the problem and for useful discussions, and is grateful to G. Chatterjee and N. Dutta for 
their support and encouragement. 
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Figure captions 

Fig. 1 : (a) A typical optimal tour is shown for = 80 cities in a unit area for dimension 
d — 2. The cities are represented by black dots and the optimal tour is indicated using a 
solid line (b) Similarly, a typical optimal tour is shown for N = 80 cities in a unit volume for 
dimension d = 3. The cities are represented by black dots and the optimal tour is indicated 
using a solid line. 

Fig. 2 : (a) Plot of the distribution p{n) of neighbours on an optimal tour oi N = 256 
cities, for d = 2. In the inset, p{n) for N = 100, in o? = 2 is plotted in the linear-log scale 
and the numerically fitted curve p{n) — 0.72 exp(—0.60n)[l + 0.05 exp(n/3.0)] is shown by a 
solid line; the error bars are due to configurational fluctuations (b) Plot of the distribution 
p{n) of neighbours on an optimal tour of TV = 256 cities, for d = 3. In the inset, p{n) 
for N — 100, in d = 3 is plotted in the linear-log scale and the numerically fitted curve 
p(n) = 0.72exp(— 0.60n)[l -|- 0.05 exp(n/3.0)] is shown by a solid line; the error bars are due 
to configurational fluctuations (c) Values of p{n) for different values of n are plotted against 
1/A^ for d — 2,to show that the frequency distribution does not have any significant finite- A" 
effect in the range 10 < A" < 256 considered. 

Fig. 3 : (a) Plot of the distribution p{n) of neighbours on a "greedy" tour of A" = 100 cities, 
ioi d — 2. In the inset, p{n) foi d — 2 is plotted in the log-log scale and the numerically 
fitted curves are shown by a dashed line (exponential decay) and a solid line (power law 
decay); the error bars are due to configurational fiuctuations (b) Plot of the distribution 
p{n) of neighbours on a "greedy" tour of N — 100 cities, for d — 3. In the inset, p{n) 
ior d — 3 is plotted in the log-log scale and the numerically fitted curves are shown by a 
dashed line (exponential decay) and a solid line (power law decay); the error bars are due 
to configurational fiuctuations. 
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Fig. 1 (a) 




Fig. 1 (b) 
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